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Abstract. Two deterministic models for Brownian motion are investigated by 
numerical simulations and kinetic theory arguments. The first model consists of a 
heavy hard disk immersed in a rarefied gas of smaller and lighter hard disks acting as 
a thermal bath. The second is the same except for the shape of the particles, which 
is now square. The basic difference in these two systems lies in the interaction: hard- 
core elastic collisions make the dynamics of the disks chaotic whereas that of squares 
is not. Remarkably, this difference does not reflect on the transport properties of 
the two systems: simulations show that the diffusion coefficients, velocity correlations 
and response functions of the heavy impurity are in agreement with kinetic theory for 
both the chaotic and non-chaotic model. The relaxation to equilibrium instead is very 
sensitive to the kind of interactions. These observations are used to think back and 
discuss some issues connected to chaos, statistical mechanics and diffusion. 
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1. Introduction 

A century after Einstein [1] and Smoluchowski [2] seminal contributions, Brownian 
motion (BM) and diffusion phenomena remain an active subject of research. Statistical 
mechanics, since its foundation, has been operating an elegant synthesis between the 
microscopic dynamical laws and the macroscopic properties of a system. In this context, 
Brownian motion is a paradigmatic example of this modus operandi. 

In the statistical mechanics framework, the minimal condition needed by 
microscopic dynamics for macroscopic diffusion can be identified in the presence of a 
mechanism leading to velocity decorrelation - memory loss. In the effort of interpreting 
BM and non-equilibrium transport in the light of modern dynamical systems theory, it 
thus comes rather natural to identify in the chaotic character of microscopic dynamics 
the main candidate to explain macroscopic transport. The instabilities of chaotic 
evolutions typically produce irregular trajectories resembling Brownian motions and 
supply a simple mechanism for memory loss. Indeed large scale diffusion has been found 
in simple low dimensional chaotic systems [H HJ [5]. This picture received theoretical 
support from the existence, in some systems, of remarkable quantitative relations 
between macroscopic transport coefficients - such as diffusivity, thermal and electrical 
conductivity -, and microscopic chaos indicators - such as the Lyapunov exponents and 
the Kolmogorov- Sinai entropy -, see e.g. [HI [71 El El [TO] . 

On the other hand, non-chaotic models generating diffusion have been proposed [HI 
[T2l [T3l [H] raising some doubts on the actual role of chaos for diffusion. However, 
these models, representing possibly the most elementary examples of diffusion without 
chaos, seem to be rather artificial: involve few degrees of freedom, often the presence of 
quenched randomness is needed together with the fine tuning of some parameters, e.g., 
for the suppression of periodic orbits [H] . Their relevance to statistical mechanics is thus 
not obvious. It is worth remarking that questions about the relevance of chaos have been 
recently raised also in the (related) context of thermal conduction problems [15], where 
non-chaotic models for heat transport were proposed and investigated [El HTJ, [THl [19] . 

It should be remarked that the above systems are non-chaotic in the sense that 
the Lyapunov spectrum is non-positive. However, there are well known examples of 
Lyapunov stable systems that display non trivial behaviours [HI [HI [161 El EE1 EH 20J. 
In the presence of dynamical randomness without the sensitivity to initial condition, 
as in quantum mechanics, an alternative definition of "chaos" or "randomness" has 
been proposed in terms of the positiveness of the Kolmogorov-Sinai entropy [21]. In 
classical systems with a finite number of degrees of freedom, as consequence of the 
Pesin's formula, the two definitions coincide. However, the proposal of Ref. [21] is an 
interesting open possibility for quantum and classical systems in the limit of infinite 
number of degrees of freedom. 

It is far from trivial to the establish the role of chaos to macroscopic transport 
properties. The problem is actually more general as chaos (in the above wider definition) 
had been often invoked to justify the whole statistical mechanics apparatus [10] and 
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Figure 1. Schematic illustration of two models used to discuss the mechanical 
diffusion: (left) hard disks and (right) hard parallel (non rotating) squares. The 
colloidal particle is displayed in red. 

to explain the irreversibility of macroscopic processes [22J. Such viewpoint coexists with 
the "more traditional" approach of Boltzmann, which stresses the role of the many 
degrees of freedom and is mathematically supported by the results of Khinchin |23j, 
Mazur and van der Linden [21] (see also Bricmont [25] and references therein). 

This paper aims to discuss the role of chaos in the context of diffusion, thus 
we compare a chaotic and a non-chaotic many degrees of freedom system, providing 
two distinct deterministic models for BM. Both models consist of rarefied gas of hard 
particles surrounding a larger and heavier particle, in the following referred to as colloidal 
particle (or colloid for brevity), impurity or test particle. As an effect of the large number 
of collisions, the deterministic motion of such an impurity behaves as a BM at long times 
(i.e. much larger than the collision time) provided its mass is much larger than that of 
the lighter gas particles. The latter condition ensures a scale separation between the gas 
and impurity dynamics. The asymptotic convergence to a BM was proved for an infinite 
one dimensional system of hard-core particles [26J and, later, in three dimensions |27j . 

The first model consists of N hard disks surrounding a larger and heavier disk 
as sketched in Fig. 1 (left). Particle interactions occur via binary elastic, hard-core 
collisions which, due to convex shape of the disk, lead to chaotic trajectories (i.e. 
exponential separation among initially close trajectories). In the following, this hard 
disk model will be denoted by the acronym HD. 

The second model, illustrated in Fig. 1 (right), consists of a gas of N hard (non- 
rotating) parallel squares (HPS) surrounding a larger and heavier square particle. In this 
case, binary hard-core elastic collisions are not central and conserve the initial parallel 
orientation of the squares. Unlike HD, the resulting dynamics is non-chaotic. This 
model is not new, the case of a HPS gas with identical particles (without the impurity) 
has been studied both in terms of molecular dynamics and kinetic theory by Frisch and 
coworkers [281 EH EH ED]. It constitutes an interesting statistical mechanics system 
characterized by the presence of an infinite number of integrals of motion, preventing 
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the system from being ergodic. Remarkably, in spite of such a pathology, the HPS gas 
gives rise to, e.g., phase transitions and transport properties as in the HD gas which, 
on the contrary, possesses ergodic and mixing properties. 

The colloid can be seen as a test particle to probe the transport properties of both 
HPS and HD. Notwithstanding the intrinsic difference in their dynamics, chaotic for 
HD and non-chaotic for HPS, their macroscopic diffusion properties are remarkably 
similar, as confirmed by the analysis of the velocity-velocity correlations, connected to 
diffusion via the Green-Kubo formula (32j [33] . Measurements of the response function 
of the velocity of the colloid and of gas particles to small perturbations confirm that the 
non-chaotic nature of HPS does not influence either the diffusion or the validity of the 
Fluctuation Dissipation Theorem (FDT). 

To discriminate the behavior of HPS from HD model one has to look at the 
relaxation to equilibrium, when the initial state is very far from it. For square particles 
relaxation properties crucially depend on the presence of the impurity, which induces 
a sort of "effective interaction" among the gas particles and allows the system to 
equilibrate to the Maxwell-Boltzmann distribution. Since collisions simply reshuffle the 
velocity components among the particles, preserving the initial velocity distribution, 
without the colloidal particle relaxation to Maxwell-Boltzmann is not possible. In spite 
of such a constraint, self-diffusion and other statistical mechanics behaviors can still be 
found [2HI ESI EH EI] (see also Sect. [3]). In the presence of the colloid, a signature of 
such a pathology survives: the relaxation time % strongly depends on the mass of the 
impurity, whilst in HD it is independent of it. Moreover, ergodicity is not fully recovered 
being the exchanges between the x and y components of the velocity forbidden. 

The paper is organized as follows. Sect. [2] introduces the two models. In Sect. El the 
numerical results for the diffusion coefficient, self-diffusion and velocity auto-correlation 
function are presented. In Sect. HI the relaxation properties of HPS and HD close-to 
and far- from equilibrium are analyzed. Sect. [5] concludes the paper with a discussion of 
the relation among chaos, statistical mechanics and diffusion based on our results. In 
the Appendix, a kinetic theory derivation of the diffusion constants for both models is 
presented in the very dilute limit. 

2. Deterministic many particle models for diffusion 

In the spirit of Smoluchovski's approach [2], it is natural to introduce a mechanical 
model for BM based on the dynamics of a macroscopically small but microscopically 
large heavy impurity colliding with many lighter particles. The large mass and size with 
respect to the gas particles allows for a separation of time scales so that the velocity of 
such an impurity is expected to follow a Langevin equation. The collisions with the gas 
provide both the friction and the stochastic kicking to the colloid |34| . Here, we focus 
on two different models in two dimensions. 
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2.1. Hard disks model (HD) 

We consider N hard disks of radius r and mass m plus an impurity consisting in another 
disk of radius R > r and mass M > m, as in Fig. 1 (left) . All particles are in a square 
box of side L, with periodic boundary conditions. Crucial parameters are the number 
density p = N/L 2 and the volume fraction ip = NTrr 2 /L 2 . In the following, we shall 
always consider very dilute systems (ip <C 1) so that the properties of the system will 
be akin to those of a rarefied gas. 

The kinetic energy coincides with the total energy H 

N 2 p2 N+1 „2 

f-,2m 2M ~ £t 2m, ' 1 ' 

and is conserved. In the second equality of Eq. (pQ) we adopt the convention that 
i = N+l indicates the colloid, i.e. Pn+i = P and m^ + i=M while irii = m for i = l, ...,N. 
Similarly, for the coordinates, we use either q { = (xi,yi) (i = 1,...,N) for the gas 
particles and Q = (X, Y) for the mass impurity, or qj for j = 1, . . . N+1 with ^at+i = Q. 
Of course, if r = R and m = M, Eq. (CE]) reduces to a system of N + 1 identical hard 
disks, which has been thoroughly investigated by accurate computer simulations [35] 
and theoretically in terms of kinetic theory of gases (see e.g. [9, 36J and references 
therein) in a variety of regimes. 

Each particle moves with constant velocity till it collides with another particle, 
event at which the velocities are updated according to the elastic collisions rule 

Pi = Pi^ ; {gij-eijjeij and Pj=Pj ■ (gy • e„)e^- (2) 

rrii + rrij J rrii + rrij 

where post-collision quantities are primed, and = Pi/mi — Pj/rrij = Vi — Vj is the 
precollisional relative velocity. The unitary vector e^, oriented as % — ► j, connects the 
centers of the two disks at contact. 

The model ([1]) simplifies for two asymptotics: for M/m — > 0, it approaches the 
Lorentz-gas model [37] , the impurity is much faster than the (now heavier) gas particles 
which can be treated as immobile obstacles; for R/r — > oo, the Rayleigh-flight model 
is recovered if the N small disks are very dilute so that most of collisions involve the 
impurity. In these two limits, kinetic theory allows for analytical approximation of the 
maximum Lyapunov exponent in terms of the system parameters |38j . 

Here, being interested in Brownian motion, we consider situations in which R^> r 
and M 3> m. In such a limit, thanks to the large time scale separation between 
the impurity and gas particles motions, it is reasonable to assume that each velocity 
component V of the colloid follows a Langevin equation (we set ks = 1) 

dV I T 

where T is the gas temperature, and rj(t) a zero mean Gaussian process with correlation 
(r)(t)r)(t')) = 8(t—t'). Of course, Eq. §S§ describes the effective dynamics of the impurity 
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for times much larger than the average collision time with the gas particles 
The friction constant 7 sets the decay of velocity- velocity correlation function 

C v (t) = (V(t)V(0)} = (V 2 )e^ = L e -* , (4 ) 

where brackets indicate time or ensemble averages. By standard kinetic-theory 
computation (see Eq. (1A.2I) where also corrections in r/R and m/M are taken into 
account) one has 



7 = 2v^^— • (5) 

Now, thanks to the Green-Kubo relation, linking the auto-correlation function to the 
diffusion coefficient D 

D = f o { V(t)V(0)}dt^ { -2 = —, (6) 
it is straightforward to derive the diffusion constant of the colloidal particle 

D c = lim-([X(t)-X(0)] 2 ) = ^- J-. (7) 

t-oo2t u w v n 1 2V27F pRl m KJ 



Notice that D c is proportional to vT and not to T as in liquids. This comes from the 
fact that in liquids the friction is temperature independent, while in rarefied gases it is 
proportional to the square root of the temperature. 

Another interesting quantity, well defined both in the presence or absence of 
the impurity, is the self-diffusion coefficient D g of a tagged gas particle. This is an 
important transport coefficient associated with the gas density field evolution. Previous 
investigations of hard disks and spheres models studied in details such a quantity as well 
as other transport coefficients (e.g., viscosity) at varying the parameters of the problem. 
In the dilute limit, the particle self-diffusion coefficient takes the form [9]: 

1 11 Ft 

D Q = lim — (\x(t) -x(0)l 2 ) = -—;= — \ — , (8) 

where x indicates the x-component of the position of a tagged gas particle. At 
high volume fractions, corrections to the formula must be considered. Although the 
expression is formally similar to ((7|), only the prefactors change, the Langevin description 
(J3J) does not apply in this case due to the absence of time scale separation. One has to 
take care to work in conditions in which the coefficient (jSJ) is insensitive to the presence 
of the impurity, so that it can be considered as a small perturbation. 

Before passing to the other model, it is important to give a warning. In two 
dimensions, as here, it is long known [391 HO] that the gas velocity auto-correlation 
function develops small, slowly decaying tails. In principle, this may lead to an ill 
defined diffusion coefficient. However, when ip — > this problem becomes relevant only 
for enormously long times [H]. Therefore, disregarding this issue is justified from the 
practical point of view and we will ignore it in the following. 
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2.2. Hard parallel squares (HPS) 

We now consider a different model of hard-core particles in which the N gas disks are 
replaced by N squares of side 2r and mass m plus a square colloid of side 2R (R > r) 
and mass M > m, contained in a box of size L x L with periodic boundary conditions. 

At the initial time, the sides of all squares are parallel (see Fig 1 (right)), and such 
parallelism is conserved (no rotation) by the elastic, hard-core collisions, which are now 
not central and amount to equal angle reflection against parallel sides. The Hamiltonian 
([1]) is still describing the system, and the collision rule <^ are replaced by 

, 2m i m j A , 2m;m j A 
Pi = Pi H ; — —{9ij ■ riijjnij and p = pj ■ — ■ n SJ n !3 , (9) 

rrii + rrij J rrii + rrij 

where gr^ is the relative velocity, and the unitary vector is directed as the normal 
to the colliding sides, that is, either along the x or y axes. At variance with the 
HD, the HPS model is not a mechanical system as its evolution does not follow the 
Newtonian dynamics, because of the constraint keeping the parallelism among the 
square sides. Since the collision rules (Q are linear, the evolution law for the tangent 
vector coincides with that of the system and the Lyapunov exponents are all equal 
to zero. Therefore, unlike the case of disks, the system is non-chaotic, although the 
presence of the corners may induce non-linear instabilities, producing a defocusing for 
non-infinitesimal displacement among two trajectories. 

In the absence of the impurity the considered system reduces to N identical parallel 
squares, which was studied by Frisch and coworkers [281 [221 [3Ql [31] for its peculiar 
properties. This system is indeed interesting under several aspects. First of all, due the 
rules of interaction ([9]), colliding particle pairs simply exchange their velocities along the 
direction of impact: similarly to Id hard rods a collision event corresponds to a relabeling 
of particles, though in the 2d case as here the relabeling is only for one component 
of the velocity. This implies that all velocity moments are conserved separately for 
the two components. It thus follows that a velocity probability distribution function 
(pdf), which is initially factorized P(v x ,v y ;t = 0) = p(v x )p(v y ), is preserved by the 
time evolution. In other words, unlike HD, the system is non-ergodic and the pdf of 
the velocities cannot relax to the Maxwell-Boltzmann distribution. However, if the 
initial distribution is not factorized, it will factorize P(v x ,v y ;t — > oo) = f(v x )g(v y ) with 
g(v x ) = / dv y P(v x ,v y ;t = 0) and g{v y ) = J dv x P(v x ,v y ;t = 0). Indeed the relabeling 
induced by the collisions decorrelate the x and y components. 

In spite of these pathologies, simulations (28] and kinetic theory computations [29] 
show that this non-ergodic system possesses many properties akin to those of the HD, 
which instead is typically considered to be ergodic and mixing. For instance, transport 
properties are well defined and in the dilute limit. The self-diffusion coefficient for an 
initial Maxwell-Boltzmann distribution can be computed [28] 
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Moreover, phase transitions similar to those found in HD can be observed for suitable 
values of the volume fraction ip = N4r 2 /L 2 [29] . 

When the impurity is present the situation changes. Due to the mass difference 
(m ^ M), many conserved quantities are destroyed and collisions can change the 
velocities pdf: relaxation to Maxwell-Boltzmann distribution is now possible. Essentially 
the impurity allows for an "effective interaction" as e.g. in [12]. However, ergodicity 
is not fully recovered: it is clear that, due to the collision rule, no mixing between 
the velocities along x and y is possible. Meaning that E x = Y,f=i~ pli/2 m i an d 
E y = Yl!i=\ v\il^ m i are separately conserved. We shall always work in isotropic 
conditions, i.e. with initial velocity pdf such that E x = E y 

In such a condition, being an equilibrium distribution well defined, we can still 
consider the limits R^> r and M 3> m and study the BM of the impurity. From kinetic 
theory (see Eq. ( 1A.5I) in the Appendix), one can compute the friction coefficient 

n-tf&ge.. (ID 

and therefore the diffusion constant 

D e =^±-M. (12) 

Comparing the above expression with (J7]) reveals that the only difference lies in a 
numerical prefactor which takes into account the different geometries of the particles. 



3. Transport Properties 

Let us now investigate, by means of numerical simulations, the transport properties 
(diffusion of the impurity and self-diffusion of tagged gas particles) in HD and HPS 
models both in the presence and absence of the impurity. 

Before discussing the results, we briefly summarize the numerical methods 
employed. For both HPS and HD, we used an event driven algorithm with the minimum 
image convention [35]. All simulations refer to the rarefied case. After several tests we 
fixed p = N/L 2 = 10 and solid fraction tp « (9(10~ 3 ), so that we can neglect all known 
corrections to the transport coefficients [36]. The diffusion constants are then measured 
in terms of the mean square displacement in the x and y direction for both the colloid 
and tagged gas particles, see Eqs. © and dHJ). Working in isotropic conditions, the 
horizontal and vertical diffusion coefficient are equal, allowing us to average over the 
two components for increasing statistics. Averages have been computed over M time 
windows in a long run. The diffusion constant D c is thus obtained by looking at 

([AQ(t)} 2 } = ^- j:[Q(t + t k )-Q(t k )} 2 ^2D c t } (13) 
■ /v fc=i 

where Q = X, Y represents any of the components of the particle position, and 
< t < (tk+i — tk)- The length of the window — t k is chosen long enough to 
observe a diffusive regime over about two decades. For the self-diffusion we proceeded 
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similarly, but working on the gas particle positions, qi = Xj,^. Being all gas particles 
equivalent, we average over them 

i Af N 

with abuse of notation the same brackets as in ({TBI are used though in (fbil) the average 
over all particles is also performed. HD and HPS simulations differ only for the collision 
rules. While computing the average displacements we also computed the velocity auto- 
correlation functions for both the colloid and gas particles Cy(t) and C v {t). The averages 
have been performed, as for the displacements, over different time windows and, for the 
gas, also over all particles. 

We work with a single impurity with R = lOr and M = 20m (tests with different 
sizes and masses have been performed). Due to the necessity of averaging over many, 
typically M = 10 3 , time windows we employ a not too large number of gas particles 
(N = 300, tests with N = 1000 have been done). Using this relatively low number 
of particles and the mass ratio M/m = 20 requires to take into account finite size 
correct iong^l Indeed the constraints of energy and momentum conservation (the latter 
ensured by the periodic boundary conditions) are known to cause a breakdown of the 
energy equipartition for simulations with non-equal mass particles [13"} I44j . 

We are now ready to discuss the results. Let us start from Fig. [2^ where we show the 
data for the diffusion and velocity correlation of the impurity and tagged gas particles 
in the case of hard disks. Although the diffusive properties of HD have been thoroughly 
studied in the literature, they serve for comparison with HPS, presented below. As 
one can see, the diffusive scaling extends over about two decades allowing for a good 
estimate of the diffusion coefficients, which are in agreement with the expectation (171) 
and (jHJ), see also the inset in Fig. [2b. The auto-correlation functions Cy and C v (shown 
in Fig. [2b) display an exponential decay. For the colloid the decay time is in agreement 
with the friction constant (j3J), making the approximation of the colloid evolution in 
terms of a Langevin equation meaningful. 

We now consider the case of N parallel squares with a mass impurity. As discussed 
in the previous section, the presence of the colloid allows for the existence of an 
equilibrium state characterized by the Maxwell-Boltzmann distribution for the particle 
velocities. However, unlike HD, the system is non-chaotic. Moreover, ergodicity is 
broken and there is no mixing between the horizontal and vertical components. As 
Fig. clearly shows the transport properties of HPS are well defined, the colloidal 
particle and tagged gas particles diffuse with coefficients as given by Eq. (f!2l and ( ITOl , 
respectively; see also inset in Fig. [3b. Moreover, the system loses memory as both the 
velocity auto-correlation function of the impurity and gas particles decay exponentially 
(Fig. [3b). At a first sight, the agreement of the self- diffusion constant with (flUj) may 

\ In practice, the temperature of the impurity will be smaller by a factor 1 — M/Mt than that of the 
gas (with Mt = J2i=i jv+i m i)- Other corrections by terms m/M and r/R should be included as well, 
see Appendix. 
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Figure 2. (a) Mean square displacement for the colloid ([AQ(t)] 2 ) (red) and tagged 
gas particles ([Aq(t)] 2 ) (blue) vs time for HD. The straight lines show the expectations 
dl|) and JHJ), respectively. The temperature is set to T = 50. (b) Velocity auto- 
correlation for colloid Cy (t) (red) and gas particles C v (t) (blue) measured in the same 
simulation as (a). The straight line has a slope given by the friction constant ([5]). Inset: 
Normalized diffusion constants RD C and rD g vs T. Solid and dotted lines correspond 
to the expectation {7} and ([5]) with the finite sample size corrections included [331 131] . 
N = 300 gas particles of radius r — 0.005 and mass m = 1 have been used, with 
number density p — 10, for the impurity R = lOr and M — 20m. Averages are 
performed over Af — 1000 time windows. Simulations on a shorter time scale with 
N = 1000 particles give the same result within statistical errors. Notice that the 
self diffusion is by definition the results of a further average over the number of gas 
particles, thus its statistics is much improved than D c implying lower statistical errors. 




(a) (b) 

Figure 3. Same of Fig. [2] for HPS. The size of the particles has been chosen so to 
have the same area of the disks. 



appear strange. Indeed, the formula derived by Frisch and collaborators [28] refers to a 
gas of identical particles, where no equilibration occurs. However, Eq. (112]) was obtained 
assuming a Gaussian distribution for the velocities, here always realized thanks to the 
presence of the impurity which allows for Maxwell-Boltzmann distribution. 

We consider now the model in which no impurity is present, previously studied 
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Figure 4. (a) (Aq 2 (t)) vs t for a gas of TV = 300 hard parallel squares with size as in 
Fig. [3] and with a Gaussian velocity distribution at temperature T = 50 and number 
density p = 10. Inset: velocity auto-correlation function. With N = 1000 we obtained 
indistinguishable results. 

by Frisch and collaborators J29J [30]. As stressed above, in such a case HD and HPS 
(i.e. chaotic and non-chaotic particle systems) are conceptually very different: while HD 
remains a well defined statistical mechanics system with relaxation to an equilibrium 
state independent from the initial conditions, HPS never reaches equilibrium, the initial 
velocity distribution remains unchanged with time. Nevertheless both systems display 
well defined transport properties. 

By computing D g for HD in the same conditions as Fig. Eh we checked whether the 
value of D g agrees with its measurement done in the presence of the colloid. We actually 
found a perfect agreement, within errors, between the two measurements (not shown). 
This confirms a posteriori that the colloid can be considered a small perturbation for 
the gas without important consequences for the transport properties of the HD gas. 
More interesting is to investigate the case of squares. 

In Fig. HI we show D g for a gas of equal hard squares with Gaussian initial 
distribution having the same temperature as in the simulations with the impurity: the 
diffusive behavior is again a robust well defined property. It is worth underlining the 
following points. First, a diffusive behavior is observed also also in the absence of 
relaxation to a statistically steady state. Second, the value of D g matches the theoretical 
value obtained by Frisch [28] and is in perfect agreement, within errorbars, with the 
equivalent quantity obtained in the presence of the impurity. Third, the velocity-velocity 
correlation function, shown in the inset, decorrelates exponentially with a good degree 
of approximation. 

The presence of diffusive behaviour in non-chaotic systems is not new. Models 
consisting of non interacting particles has been already considered in Refs. [T2l [13]. 
However, such models, though interesting from a dynamical system point of view, 
consist of independent particles and thus are somehow far from a statistical mechanics 
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perspective. The HPS system here investigated is non-chaotic but having many degrees 
of freedom in interaction it constitutes a valuable statistical mechanical system. 

We conclude this section by stressing that, although the presence of the colloid 
changes conceptually the properties of the system, transport properties in the presence or 
absence of the impurity are quantitatively the same, provided the Gaussian distribution 
is chosen. This means that, at least, for transport properties the colloid can be 
considered as a small perturbation to the system, but for a non generic state. 



4. Relaxation Properties 

The analysis of the relaxation processes associated with spontaneous or induced 
statistical fluctuations represents the classical approach to probe the macrostates 
explored by a system. For instance, Fluctuation Dissipation Theorems (FDT) |33j, 
relating the behavior of spontaneous fluctuations at equilibrium to the average response 
of a system to infinitesimal perturbations, establish a connection between equilibrium 
(correlation functions) and non-equilibrium (response functions) quantities. We carried 
out a set of simulations to determine the relaxation properties of the system both close-to 
the equilibrium state, characterized by a Maxwell-Boltzmann distribution, and far-from 
it (e.g., starting from a uniform distribution) to probe the possible influence of chaos 
on the statistical properties. 



4-1. Close to equilibrium 

To gain information about transport properties by studying the relaxation to equilibrium 
it is useful to introduce the response function to small impulsive perturbations applied 
to a component V of the velocity of the impurity. This can be obtained, e.g., by applying 
a force fit) = FS(t) which acts only at t = with F < 1. The result of fit) is to 
cause an instantaneous (very small) variation of the velocity V(0) — > V(0) + 5Vo, with 
SVq = F/M. One can thus define the average response function as Rv(t) = (SV(t)) e / '5Vq, 
where ([. . .]) e denotes an ensemble average at fixed time in the presence of an impulsive 
perturbation. The response Ryif) bears important information on the transport 
properties of the system being. Indeed, if the velocity distribution is Gaussian the 
classical FDT relation [33] tells us that Rv{t) coincides with the normalized velocity 
correlation function 

= as) 

This can be seen as the differential form of the Einstein relation connecting the 
asymptotic speed of the particle to the mobility under the effect of an infinitesimally 
small force. In the following, we present the measurement of the response function 
Ry{t), which will be compared with the correlation Cy(t) to test whether FDT holds. 
Similarly, we analyze also the response function of a tagged gas particle Ry{t) when the 
impurity is absent. 
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Figure 5. (a) Rv(t) vs t measured for SVq = a^jT/M with a = 0.02, 0.05, 0.1 as in the 
legend and CV(f)/CV(0) (solid line) computed in the same conditions. The parameters 
are as in Fig. [2] with T = 50. Note the good collapse, within statistical errors, of the 
response functions for a = 0.02 and a = 0.05 and the correlation function. Deviations 
can be appreciated for a = 0.1 (presumably for such a value of the perturbation one 
exits the linear response regime), (b) the same of (a) for HPS system with parameters 
as in Fig. [3l the agreement between the decay of the correlation and the response 
function is also in this case very good. For both (a) and (b), the three curves have 
different length because decreasing the initial perturbation value the signal is spoiled 
by noise at earlier times. 



In order to numerically compute the response function, we adopted the following 
protocol. Consider a system HPS or HD in the presence of the colloidal particle and let it 
evolve till the equilibrium state is reached. At this time, that we call t = 0, the velocity 
of the colloidal particle is perturbed by a small amount V(0) — > V(0) +5V . In principle, 
the perturbation should be infinitesimal (i.e. SVq — > 0). This is however infeasible in 
practical computations, we then considered three different perturbation values defined 



in fraction of the root mean square velocity SVq = a^jT/M with a = 0.02,0.05 and 
a = 0.1. In this way, comparing the different numerical experiments, we can test 
a posteriori the validity of the linear response theory. This small perturbation on the 
velocity is expected to be re-adsorbed, meaning that the velocity of the impurity should, 
after a while, assume values drawn from the equilibrium distribution. The procedure is 
thus repeated for many (typically 10 4 — 10 5 ) times. The time history of the colloidal 
particle V(t) is followed in each experiment so to obtain its average evolution (V(t)) e , 
from which the response function Ry(t) can be defined as 

p ,,x (V{t))e hR x 

dVo 

note that we used that (V(0)) e = 0. 

Equation ffl6]) has been verified in our simulations for both HPS and HD, and the 
comparison between the numerical results is shown in Fig. [51 One should notice first that 
for both HD and HPS the response measured for different perturbations superimpose, 
meaning that we are in the linear response regime. Moreover, the fair superposition 
of the exponential decays of Cy{t) and Rv(t), up to statistical errors, constitutes the 



Transport properties in chaotic and non- chaotic many particles systems 



14 




Figure 6. The same of Fig. EJa) showing the FDT relation for the HD (a) and HPS 
(b) gas particle (in the absence of the impurity). For HPS the initial distribution of 
the gas velocities was chosen Gaussian with temperature T = 50 so to compare with 
the correlation function as measured in Fig. [4] 



numerical evidence for both systems to obey FDT relation. 

The conclusion that can be drawn is that whenever the chaotic HD and the non- 
chaotic HPS are prepared into an equilibrium state compatible with the thermodynamic 
parameters T and p, the behavior of Rv{t) does not reveal any difference between the 
two systems. 

The above procedure can be applied to a tagged particle so to define the response 
for gas particles R v (t) which is connected to the correlation function C v (t). We perform 
such a measurement for the system of HD and HPS without the impurity. Being the 
perturbation very small the measurement is still meaningful though the system is not 
equilibrating. 

Figure [6] shows the average response functions for a tagged particle in both HD 
and HPS, together with the comparison with the correlation functions. There is a fair 
agreement between R v (t) and C v (t) / C v (0) for both HD and HPS models. Of course, in 
the case of the HPS the origin of the validity of a FDT relation cannot be ascribed to 
the presence of chaos. Clearly it can only come from the presence of many degrees of 
freedom and should have a probabilistic origin. 

4-2. Far- from equilibrium 

We consider now the case of relaxation when the system is prepared into a state far-form 
equilibrium, for which the difference between the two models becomes evident. This 
can be understood from the outset by recalling that HPS systems with identical squares 
cannot relax due to the collision rules that merely relabel the velocity components. 
Unlike HPS, hard disks collisions, also thanks to chaos, mix the velocity components at 
each impact and allow for a fast relaxation to a Maxwell-Boltzmann distribution. 

When the colloid is introduced, relaxation to equilibrium becomes possible also in 
HPS because impacts against the impurity break the relabeling process and provide a 
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mechanism for the transfer of energy (at least, separately for the x or y components 
which do not mix), this is similar to the problem of the adiabatic piston considered in 
Ref. [32] • However, since only collisions with the impurity contribute to the relaxation 
process of the gas, the time-scale for reaching the equilibrium state crucially depends 
on its mass M and size R. This contrasts with the HD model, for which the time scale 
of relaxation is essentially unaffected by the characteristics and/or the presence of the 
colloid. 

In our simulations, we prepare the HD or HPS gas in a state with a spatially 
homogeneous distribution of particles and flat velocity distribution, i.e. P(v x ,v y ) = 
g{v x )g{vy) with g(v) = l/(2t> ) for \v\ < vq and zero elsewhere. The value of vq is 
fixed by imposing the temperature T of the system, i.e. Vq/6 = J dvg(v)v 2 = T/m. The 
colloidal particle is initialized with random velocity extracted from the same distribution. 
The system is then let evolve under the event driven dynamics, and the velocity pdf 
monitored by computing 

iv x (tf 3v,(tf 

The symbol [. . .] indicates the average at a given time t over the gas particles, i.e. 
v% = l/NJ2iLi v 1i- At t = 0, JC = K = 3/5 while for t — > oo the Gaussian result 
/C = = 1 should hold, being the system relaxed. To have a smooth behavior we 
average K(t) over many independent runs. We thus obtain ()C(t)) e which is shown in 
Fig. [T^i for HD and is well described by the fitting function 

(K(t)) e = K oc + (K - ^ 00 )e-' /r " (17) 

where the fitting parameter %. provides an estimate of the relaxation time to equilibrium 
of the system. The brackets ([. . .]) e denote averages over the realizations. 

The perfect collapse of the curves (JC(t)) e for the HD system (Fig. [7h) obtained for 
four values of the ratio M including the case M — m — 1 (i.e. without the impurity), 
clearly indicates that the relaxation process of the disks is independent of the colloid. 
Unlike HD, for HPS the mass M of the impurity is crucial in determining the relaxation. 
This is shown in Fig. [7b, where we report the behaviour of fitted by using (fTTI) . as a 
function of M. As one can see, % diverges for M — > m (identical squares) and M — * oo 
(immobile impurity which again does not allow for the exchange of energy, leading to 
the impossibility of relaxation). Notice that asymptotically % seems to grow linearly 
with M. 

This results, although obvious when considering the different nature of the 
collisional processes occurring in the two systems, confirm that the relaxation to the 
equilibrium of HPS is much slower than that for HD and crucially depends on the 
impurity characteristics. 
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Figure 7. (a) (IC(t)) e vs time for HD with different impurity mass values including 
M = m (i.e. absence of the colloid). All curves collapse confirming that in the case 
of disks the relaxation time is independent of the presence of the colloid. The black 
curve represents the fit (fT7| . (b) Relaxation time %, vs M needed for HPS to relax 
from a uniform to the Gaussian equilibrium distribution at T — 50. Other parameters 
are fixed as in Fig. [3l The dotted straight line suggests a compatibility with a linear 
behavior at large mass ratios. The red horizontal line indicates T r for the HD, which 
is put here for comparison and to show the independence of relaxation process on the 
presence of the colloid as already evidenced in (a). 



5. Final remarks 

From the results presented in this paper we obtain good indications that a non-chaotic 
system with many degrees of freedom, consisting of a gas of hard squares, provides a 
suitable model for Brownian motion which is equivalent, at least at a simulation level, 
to the corresponding chaotic model, where squares are replaced by disks. 

A deep understanding the role of chaos for Brownian motion and for transport 
properties is a very important issue deserving further comments. With reference to 
previous works, we now discuss this issue by stressing how the presence/absence of 
ergodicity, mixing and chaos in the microscopic dynamics (may) influence the statistical 
mechanics of macroscopic systems. 

Let us start with a few remarks on equilibrium statistical mechanics, where the 
problem of connecting (micro) dynamics and (macro) statistical features is at the origin 
of Boltzmann's ergodic hypothesis. First, it is worth noticing that ergodicity, in its 
strict mathematical formulation, is an extremely demanding property. Second, in spite 
of its theoretical importance, it is not completely satisfactory from a physical point of 
view as it involves global asymptotic limits, rarely encountered in practice. However, for 
the foundation of statistical mechanics, a widespread consensus exists on the key role 
played by the huge number of degrees of freedom involved in a macroscopic system rather 
than ergodicity. This point of view received mathematical support from the works by 
Khinchin [23], Mazur and van der Linden [23], and others. They proved that statistical 
mechanics works independently of ergodicity thanks to the existence of meaningful 
physical observables (the so-called sum functions) which are nearly constant on the 
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energy surface, apart from regions of vanishing measure. Support to this picture comes 
from the results of Frisch and coworkers [28, 29, [30], [31], who found robust statistical 
phenomena in a trivially non-ergodic system. Nevertheless, Khinchin's statements 
cannot be the ultimate grounding of statistical mechanics because not all physically 
important observables belong to the class of the sum functions. 

Since chaos grants the validity of some "statistical laws" even in few degrees of 
freedom systems, one could be tempted to invoke it as the sufficient ingredient to build 
a robust statistical mechanical approach grounded on Hamiltonian systems. However, 
for this to be true, chaos in the microscopic dynamics should be enough "strong". 
Indeed, results^ of extended simulations in high dimensional systems [I6l Wf[ HE] have 
shown that chaos may be not enough to ensure the validity of the equilibrium statistical 
mechanics. 

Let us now discuss this issue in the non-equilibrium statistical mechanics context, 
where the analogous of ergodicity is the mixing condition. As before, this condition 
is very demanding, being related to the T-space (the set of positions and momenta of 
system particles); while the study of macroscopic systems usually focuses on physical 
observable involving some projection procedures, amounting to neglect (or average) the 
effects of a large number of degrees of freedom in favor of few relevant variables. For 
instance, in the case of elementary transport properties the observable that matters refer 
to single particle properties: mean square particle displacement, correlation function and 
the response function either of the colloidal particle or the single gas particle. Therefore, 
one can wonder about the microscopic conditions ensuring a "good" statistical behavior 
for the above quantities. 

We can start by considering few degrees of freedom systems, where also simple 
deterministic chaotic models may exhibit transport properties similar to those of more 
realistic systems. Paradigmatic examples are chaotic billiards and the Lorentz gas, 
where particle trajectories are chaotic as a consequence of the convexity of the obstacles. 
Numerical and theoretical works have shown that, in these systems (under appropriate 
hypothesis such as hyperbolicity etc.) the transport coefficients can be quantitatively 
related to chaos indicators [HI [7]. This would suggest that chaos is tightly related to 
transport. However, several examples of non-chaotic deterministic systems, such as a 
bouncing particle in a two-dimensional billiard with polygonal randomly distributed 
obstacles, possess robust transport properties [121 H31 EEU EES1 [EE] • For these non-chaotic 
models, it has been proposed that a sort of non-linear instability mechanism is required 
to observe diffusion [T3]. The existence of non-chaotic models able to display diffusion 
poses some doubts on the possibility to make strong statements on the role of chaos for 
transport. It is however important to stress that in all these models the particles do 
not interact, and therefore, at least from a statistical mechanics point of view, they are 

§ In chaotic high dimensional systems one may have a sort of "localized chaos" without a globally 
irregular dynamics. This phenomenon is somehow the high dimensional analogous to the presence 
of chaos in bounded regions with the absence of large scale diffusion observed in low dimensional 
symplectic systems in situations below the resonance overlap [45j . 
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rather artificial. 

More interesting it is thus to consider many degrees of freedom systems, such a 
those investigated in this work. In this case the existence of quantitative relationships 
among chaos indicators and transport coefficients is, to the best we know, less clear 
(see for example [IHl [50]). Nevertheless, chaos has been proved to be relevant in the 
establishment of some non-equilibrium properties [51] . Moreover, it is fair to say that 
diffusion of particles is rather common in chaotic many body particle systems. On 
the other hand, the non-chaotic model here investigated together with the previous 
results by Frisch and coworkers J28J I2H1 SO] indicate that transport properties for 
both impurity and gas particles agree with the prediction of kinetic theory and are 
indistinguishable from those of the (mixing) hard disk model. Therefore chaos, at least 
in the sense of positive Lyapunov exponents, cannot be invoked to explain the observed 
statistical behaviors. Nevertheless, as for low dimensional models [H], also in HPS a 
non infinitesimal mechanism of instability can be induced by the presence of singular 
corners of the squares, and these likely play a role for the diffusive behaviour. 

We interpret these findings in the framework developed by Khinchin: the "good 
transport" properties observed in the non-mixing system result from the large number 
of particles and not from chaos. This is well evident for the correlation and response 
functions for the hard squares when, e.g., all particle are identical and the collision 
dynamics reduces to a mere relabelling. In such a case the exponential relaxation 
of C v (t) and R v (t) is just a probabilistic consequence of the exponential distribution 
of the time interval between two consecutive collisions. We stress that here diffusive 
properties are the outcome of the action of many degrees of freedom, and not of non- 
linear instability mechanisms as in low dimensional chaotic and non-chaotic models. Of 
course, as discussed in Sect. I4.2[ chaos may favour the equilibration of the system. 

As a last remark we note that the Gallavotti-Cohen [52] fluctuation theorem seems 
to apply to non-chaotic models, at least in finite time intervals, as shown by Benettin 
et al. [53] who investigated a non-equilibrium version of the Ehrenfest wind-tree model, 
which is non-chaotic. In such a system, although the maximum Lyapunov exponent is 
zero, the presence of long irregular transients, introduces an "effective randomness". 

In conclusion, we think that it is very difficult to decipher the signature of chaos 
in transport phenomena observed in many particles systems because, as shown in this 
paper, it can be overwhelmed by the emergence of an "effective dynamical randomness" 
due to the combination of: a) coarse-graining procedure, b) finite scale instability and 
c) presence of a huge number of degrees of freedom. The characterization of this 
"effective randomness" requires the renounce to asymptotic limits (arbitrarily long time 
and arbitrary resolution) in favour of a finite time and/or finite resolution analysis 

[mini Ea E5]. 
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Appendix A. Computation of the diffusion coefficient 

In this Appendix, we detail the computation of the diffusion coefficient for the colloidal 
particle in an uniform and rarefied HD and HPS gas, following elementary kinetic 
theory. The basic idea is to estimate the average drag force exerted by the gas particles 
which collide with the impurity, by calculating the average exchanged momentum in the 
collisions. 

Consider a rarefied HD gas at equilibrium, and focus on the collision of the colloidal 
disk characterized by its mass M, radius R and precollisional velocity V with the gas 
particles which are characterized by m,r,v, respectively. According to Eq. (j2J), the 
implulse transferred in the collision is 



g = v — V being the precollisional relative velocity. The rate of such collisions can be 
obtained by considering the equivalent problem of a colloid, at rest, with radius r + R, 
and hit by a flux of pointlike particles moving at relative velocity g. The rate is then 
determined by counting the number of pointlike particles hitting the unit surface per 
unit time for a given orientation e. This number corresponds to the particles contained 
in the collisional cylinder of infinitesimal base (R + r)dd and height p\g ■ e|0(— e ■ g)5t, 
as shown in Fig. IA11 The unitary step function 0(s) selects the condition, g ■ e < 0, to 
have a collision. 

Accordingly, the mean impulsive force in the normal direction e = {cos(#), sin(#)} 
selected by the #-angle that e forms with vector V (taken as x-axis direction) 



the average is meant over the equilibrium distribution of the gas velocities P(v) = 
m/{2nT) exp[— m\v\ 2 /(2T)]. It is convenient to make the change of variable (v x ,v y ) — > 
(g x ,g y ) and then perform the approximation P(v) ~ P(g)[l — m/Tg ■ V], justified in 
the limit M ^> m. After this manipulation, the integral (lA.lj) becomes 



MAV = M(V' - V) 



2mM 



m + M 




2mM 
m + M 



p(R+r) JJdO JdvP(v)Q(-e-g)\g-e\(g-e)e, (A.l) 



( 




4m 2 M 



p(R + r)J o W d6 JdgP(g)Q(-e ■ g)(g ■ ef{g ■ V) e 



T(m + M) 



This expression, recast in the form MV = — 7MV, allows to explicitate the friction 
coefficient 7. Passing to polar coordinates, g x = gcos(a),g y = ^sin(a), simplifies the 




Figure Al. Cartoon of geometric construction used to compute the collision rate 
between an HD and HPS colloidal particle characterized by (M, R, V), and gas particles 
(to, r, v). The equivalent problem considers the colloid at rest, but with increased size 
R + r, and pointlike gas particles moving at the relative velocity v — V. The collision 
rate amount counting the number of gas particles contained in the collisional cylinder. 

integral (a being the angle between g and V, whose direction coincides with x-axis). 
Indicating by <ft the angle between e and g, we end up with 

/■27T /"OO rl'K 

I d6 I dgP(g)g 4 / da\ cos(0)|cos(0)6[— cos((f))]V cos(a) cos(6 ) ) . 
Jo Jo Jo 

where the angles a, 6, <p are related by 9 — a = n — <fi. The integration over g yields 

3/ (27r)y / 7r/2(m/T) -3 / 2 , and that one on the angles gives the value 47r/3. We then obtain 

the friction coefficient 

7 = 2 ^gg^f 1+r/ M. (A.2) 



M \l + m/M 

where the last factor takes into account finite size and mass corrections. By using Eq. ([6]) 
the diffusion constant is obtained as 

T 11 I? (l+m/M\ 

~ M7 _ U^^R^m \ 1+r/R J (A ' 3) 
The same computation can be repeated for HPS, the only difference lying in the 
geometry of the problem (Fig. IAII) . Thanks to the symmetry, it is convenient to 
decompose the average impulse transferred from the gas to the colloidal particle in 
the x and y components (along the directions (1,0), or (0,1)). For example, in the 
x-direction we have 

8V \ 4m 2 M r 

M ~bt ) ~ ~ T(m + M) p{R + r)Vx ^ ± J d ^P{9 x m-9 x n x ){g x n x ) 2 \g x \ 

Notice that the sum over n x = ±1 stems from the two identical contributions to the 
exchanged momentum given by the left and right collisions occurring on the opposite 
sides of the square. The integral can be solved in Cartesian coordinates, and considering 
that the constraint imposed by 0[— g x n x ] amounts to a factor 2. We thus obtain 

, r SV x \ 8m 2 M .„ , T , r°° , „. , , 

M ir) = _ T(m + M) + r)V - L d * M * ' 
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The integration gives the value v / 2/^(m/T) 3 / 2 yielding the final results for the friction 
coefficient 

[2 Vm P RVT f l + r/R \ 
7 Vvr M {l + m/Mj' 1 ' 

and thus the diffusion coefficient reads 

n _ A 1 [T( l+m/M \ 
c ~ 8V2 P R\m{ 1 + r/R ) ' ( 5j 
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